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■ Abstract 

o\ : 

Q\ , One-loop quantum corrections to the classical vortices in 2+1 dimensional 

I 0(3)-models are evaluated. Skyrme and Zeeman potential terms are used 

■ to stabilize the size of topological solitons. Contributions from zero modes, 
>—i . bound-states and scattering phase-shifts are calculated for vortices with wind- 

i ing index n = 1 and n = 2. For both cases the S-matrix shows a pronounced 

■ series of resonances for magnon-vortex scattering in analogy to the well- 
. established baryon resonances in hadron physics, while vortices with n > 2 are 

J> I already classically unstable against decay. The quantum corrections destabi- 

' lize the classically bound n = 2 configuration. Approximate independence 

■ of the results with respect to changes in the renormalization scale is demon- 
ic i strated. 

o . 
o\ ■ 

^ ; PACS 12.39.Dc, 12.39.Fe, 75.30.Ds, 75.50.Ee 

p Keywords: skyrmions, 0(3)-model 

^ ; 1 Introduction 

I Effective field theories have found increasing interest as powerful tools for describ- 

^ ■ ing the dynamics of physical systems where global symmetries are spontaneously 

. broken and continuous order parameter fields represent the relevant low-energy de- 

grees of freedom. Depending on space dimensionality and the manifold on which 
the fields live classical localized static solutions may fall into topologically distinct 
classes characterized by integer winding numbers, carrying energy, momentum and 
internal properties which suggest their interpretation as particle-like excitations of 
the uniform ground state. Their spatial size is determined by scaling properties of 
different competing terms in the effective lagrangian. 

Quantization of effective field theories allows to assign proper quantum numbers 
to quasi-particle properties, identify excited states of quasi-particles, and calculate 
loop corrections to the classical results for observable quantities, which may be cru- 
cial for experimental verification. The evaluation of loop corrections necessarily 
brings about the need for renormalization. Although, generally effective field theo- 
ries might be non-renormalizable in the strict sense, they still may be renormalized 
order by order in terms of a gradient expansion. A well-known example is chiral 
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perturbation theory (ChPT), applied successfully in hadron physics. In the vacuum 
sector such an expansion is truncated by allowing only for external momenta small 
compared to the underlying scale of the theory. However in the soliton sector, the 
soliton itself constitutes "external" fields with gradients comparable to the scale of 
the theory which cannot be made small by definition. Therefore, in dealing with 
solitons, the problem of truncating the expansion can only be solved by the ad hoc 
assumption that the renormalized couplings of the higher gradient terms are small. 

A prominent example for the application of this program are 3D-SU{Nf) skyr- 
mions in A^/-fiavor meson fields, where the topological charge is identified with 
baryon number, with impressive success for baryonic properties, resonances, and 
meson-baryon dynamics. Similarly, the conjecture to consider 2D-0(3) spin textures 
as charged quasi-particles in ferromagnetic quantum Hall systems and antifer- 
romagnetic high-temperature superconductors , with topological winding density 
identified with the deviation of the electron density from its uniform background 
value, suggests a corresponding investigation. For magnons and vortices with unity 
charge such attempts have been presented for ferromagnets and antiferromag- 
nets 0] . The field theoretical approach forwarded here is rather related to the anti- 
ferromagnet due to its preserved time-reversal invariance which makes the analysis 
fairly similar to what applies to relativistic systems. The ferromagnet where time- 
reversal invariance is broken would require the consideration of the Landau-Lifshitz 
dynamics. The evaluation of the Casimir energies involves bound-state energies 
and a sum over scattering phase-shifts [^. This provides also complete informa- 
tion about resonant excited states in the continuum which in the 3D-SU{Nf) case 
successfully describe well-established baryon resonances. 

For the general outline of the necessary steps and technique we discuss the case 
where the second-order exchange energy is taken in the time-reversal invariant form 
of the non-linear sigma model as it naturally appears in relativistic field theories. In 
two dimensions it is scale invariant and therefore irrelevant for the spatial extent of 
the static localized solution. To fix the soliton size two more terms are needed. We 
use the standard fourth-order Skyrme term, and a symmetry-breaking coupling of 
the 0(3)-order parameter to an anisotropy field. The static part of the Skyrme term 
represents a local approximation to the Coulomb energy of the charged excitations; 
in our field-theoretic example we keep also the time-derivative parts of it. So our 
discussion will mainly serve a demonstrative purpose as a model field theory. 

We shall specifically address the question of the binding energy of doubly charged 
quasi particles. For 3D-SU (Nf) skyrmions this question has a long history since it 
was discovered that in the winding number B = 2 sector the static lowest-energy 
solution is bound with respect to decay into two separate B = 1 skyrmions and dis- 
plays only axial symmetry in contrast to the radially symmetric hedgehog skyrmions. 
There has been much discussion about the physical relevance (in the deuteron) of 
such a torus configuration, and it was speculated that quantum corrections might 
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reverse the sign of the binding energy E{B = 2) — 2E{B = 1). Evaluation of loop 
corrections to the 3D-SU{2) B = 2 torus is a formidable task which to our know- 
ledge has not yet been achieved. It is interesting that the same situation occurs for 
the 2D-(9(3) skyrmions: the classical n = 2 solution shows a ring-like density distri- 
bution and it is bound with respect to decay into two individual B = 1 skyrmions. 
In this case, however, the evaluation of loop corrections for both configurations is 
of comparable complexity, and we will show that they in fact reverse the sign of the 
binding energy. 

Actually, there is an ongoing discussion P, |^ about the binding of 2-skyrmions 
also in condensed matter applications. However, it should be stressed that for 
ferromagnetic systems the time- dependent part of the effective lagrangian has to be 
replaced by the T-violating Landau-Lifshitz form with only one time derivative. For 
antiferromagnets the time derivative coupling of the staggered spin to the magnetic 
field, which does not contribute to the static stabilization should be included |^, 



10]. And, in any case, for quantitative conclusions, the non-local character of the 



Coulomb energy should be respected [O 



2 Static solutions 

The Lagrangian of the 0(3) model in 2 + 1 dimensions is conveniently written in 
terms of the 3-component field $ satisfying the constraint = 



The three terms represent the non-linear sigma (Nia) model, the Skyrme and the 
potential term. There exists a conserved topological current 

T^' = —et^^P^^d^^ X dp^) , 9^r^ = . (2) 
Stt 

The radially symmetric hedgehog ansatz, written in polar coordinates (r, if ) 



f sin F{r) cos nip 

sin F(r) sinnip I , (3) 



\ cosF(r) 
corresponds to winding number 

/ rfVT° = -[cos F(oo) - cos F(0)] =n, = (4) 

J 2 Aur 



and the soliton's size may be defined according to 

- fd'rr'T' = - — 
n J 4:Tr 



<r^ >^=- f d^r r'^T^ = I d^r rF' sin F . (5) 
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It is convenient to absorb the length Xj^fjem which sets the scale for the size of 
localized structures into the dimensionless spatial coordinate x = y/Jemr. The 
static energy functional connected with the lagrangian (|l]) is then given by 



e: 



class 
n 



P 

2 



d X 



x^ 



+ a — + 2a 1 - c 



(6) 



with the abbreviations s = sinF and c = cosF. Technically, a = m/fe is the only 
nontrivial parameter, while 47r/^ sets the overall energy scale. We shall therefore 
present results for different values of a. The limit a —>■ describes the pure Belavin- 
Polyakov (BP) solution |]12|. The opposite limit, a —>■ oo, is technically also well 
defined; it describes a system without spin-spin aligning force. 

The Euler-Lagrange equation following from the variation of the energy func- 
tional 
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(7) 



is a second order non-linear differential equation which is solved subject to the 
boundary conditions 

F(0) = 7r, F(x-^oo) = 0. (8) 

In principle the boundary condition at the origin according to (^) could be F{0) = 
with u = 1,3,... an odd multiple of vr, however it turns out that the solution with 
u = 1 has the smallest static energy. 



PROFILE 

a=0.316 




ENERGY DENSITY 

a=0.316 




Figure 1: Hedgehog profiles and energy densities for the skyrmions with n = 1 and 
n = 2. The parameter = 0.1 was used. 



The hedgehog profiles and the corresponding energy densities for = 0.1 are 
shown in Fig. 1 for the 1- and 2-skyrmions. It is noticed that the maximum energy 
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Table 1: Classical soliton energies and square radii for 
the stable 1- and 2-skyrmions for various parameters a. 



a 


0.0 0.01 0.0316 0.1 0.316 1.0 


El 
El 


1.00 1.03 1.08 1.21 1.56 2.57 
2.00 2.03 2.10 2.30 2.94 4.85 


< >i 

< >2 


oo 2.44 2.11 1.80 1.57 1.42 
TT 3.11 3.06 2.99 2.78 2.76 



density for the n = 2 soliton is not located at the origin but in a spherical shell of 
radius a; ~ 1. 

The classical energies and square radii are given in Table [^. In the two limiting 
cases, a ^ and a ^ oo, the differential equation (|^) may be solved analytically 
(see appendix). The corresponding static energies and radii are 



E^ = n, < >n= a/^M - 1) ■ a«l 

™ ' V 3 ^ ' sm(7r/n) 

E° = 4n/3-a, < >„= 4n/3 a»l. (9) 

It is noticed that for the BP solution with a — the 1-skyrmion's radius diverges 
which is due to the particular choice of the potential term in (|l]). This divergence 
is perceptible only for extremely small values of a, for a > 0.001 the radius for the 
1-soliton is still smaller than that of the 2-soliton which stays finite in the limit 
a — *• (for a more thorough discussion of this problem we refer to the appendix). 
The dependence of the energies on the dimensionless parameter a is shown in Fig. 
2. Independently from this parameter E2 < 2Ei always holds, the 2-skyrmion 
is classically stable against decay into two 1-skyrmions. Hedgehog solutions with 
higher topological charge, n > 3, are not stable (see subsection 3.3 for their stability) 
and do not represent the minimum energy configuration of the corresponding sector. 
The 3-soliton's minimum energy configuration is not rotationally symmetric, it is a 
linear molecule which consists of three distorted 1-solitons. Higher multi-solitons 
are then supposedly molecules of stable 2-skyrmions (and 1-skyrmions in case of n 
odd) |]l3l. This picture may change fundamentally when the Casimir energies are 
taken into account. 

3 Small amplitude fluctuations 

The evaluation of quantum corrections is based on the normal modes of the classical 
configurations. With the appropriate boundary conditions these normal modes de- 
scribe excited locahzed (bound) states, and the scattering of the vacuum fluctuations 
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Figure 2: Classical soliton energies for n = 1 and n = 2 as functions of the 
dimensionless parameter a. For better comparison the energies plotted are divided 
by the corresponding topological charge. The 2-skyrmion is classically stable against 
decay into two l-skjnrmions. 

(mesons or magnons) off the soliton backgound. (The S-matrix for magnon-vortex 
scattering has been considered (without the stabilizing Skyrme and potential terms) 
in 



14 1 for antiferromagnets and in [|T5| for ferromagnets). Because of the constraint 
$ = 1 the model effectively possesses two independent fields. We therefore 



introduce two-component small amplitude fluctuations 



ri = VnVL + <PnVT : 



COS nip 
sin Tvp 



— sm nip 
cos mp 



(10) 



which we decompose into components parallel and perpendicular to the soliton con- 
figuration with winding number n. Of course, there are many different ways to 
parametrize the total time-dependent field $ which lead to different equations of 
motion for the corresponding fluctuations. However, the bound-states and the 
scattering matrix, in particular the phase-shifts are independent of the chosen 
parametrization. A very convenient choice different from the common Polyakov 
parametrization [ff6l, which only at first sight may seem complicated, is given by 



f„ sinF(l - 77V2) + cos Fr]L + <Pj1t 
cosF(l — rf/2) — sin Fr]L 




(11) 



The main advantage of this parametrization is that it leads to a fiat metric for 
the non-linear sigma model part of the lagrangian. For 3D-SU{N) skyrmions this 
parametrization is well-known as Callan-Klebanov ansatz ITi 
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It is now straightforward to write down the e.o.m. for the fluctuating components 
TjL and TjT- The Lagrangian (|l|) has to be expanded to quadratic order in the 
fluctuations, then the e.o.m. can be read off. We give these coupled linear differential 
equations explicitely, where we have already exploited the time-dependence as well 
as the angular dependence of the fluctuations 

r/L = E /M(x)e^*^'^e-^-* , m = 9m{x)^^'' e-^^^' • (12) 

M M 

The e.o.m. decouple in the magnetic quantum number M, which is the analogon 
of the phonon- (or grand-) spin familiar from scattering calculations for 3D-S'[/(A^) 
skyrmions 

^ , t. ,1 ^< o2^^ 2n^c(F's\' 
— {xolIpi) H TJM + ^[ot- 2s )/m - a /m + acfM 

X X^ X^ X \ X 



nMc.^ nMF's , 2nM f F'sV 



1 + bT)gM + a — + a gn = hif. 



9 V ' -'^/.-?-'*'^' O 



M 



I, , M2, ^ .„ n^c^ r?c(F's\' , , 

— \^9m) + —OTgu -{F —)gM - a gu + acgu (13) 

Ou Ob Ob Ob \ OC I 

nMc.^ , , „ nMF's., nM ( F' s\ .. 
^-(1 + OT)tM - a — fM + a Jm = uj brgM ■ 

X'^ X \ X j 

Here we introduced the metric functions 6l(x) = 1 + an'^s^ /x^ and bT^x) = 1 + aF''^. 
The energies u are understood in natural units of y/Jem such that the threshold 
occurs at u; = y/a. The equations are invariant with respect to the simultaneous 
replacements M —>■ —M and gM —>■ —gn- These equations contain all the informa- 
tion about zero-modes, bound-states, the scattering matrix and the stability of the 
hedgehog solutions. 

3.1 Zero- modes 

The soliton's energy (|]) is invariant under spatial rotations around the 2;-axis (which 
for the hedgehog is equivalent to an iso-rotation around the internal 3-axis) as well 
as under a translation in the x-y plane. The infinitesimal rotation (iso-rotation) 
and translation give rise to zero-modes which are solutions of the e.o.m. (0) for 
^2 = 0. 

The rotational zero-mode is obtained by a shift of the azimuthal angle ^ (p + a 
in the hedgehog ansatz (^. Comparison with (|1^) and (|1^) shows then that this 
zero-mode 

/o(x)=0, g^i^x) = s (rotational zero-mode) (14) 
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is located in the M = partial-wave. Similarly, the translational zero-mode r 
r + a 

ns 

fi{x) = F' , gi{x) = — (translational zero-mode) (15) 

sits in the M = 1 partial-wave. 

With the stability condition (|^) it is easily checked that the above modes satisfy 
the e.o.m. for a;^ = 0. As both zero-modes possess a finite norm, they will act as 
bound-states in the scattering calculation and influence the phase-shifts via Levin- 
son's theorem. Besides these zero energy bound-states there are also bound-states 
at finite energy. 

3.2 Bound-states 

In both systems, n = 1 and n = 2, we observe bound-states in the M = n partial- 
wave. With decreasing parameter a the energies of these bound-states move towards 
the threshold, in particular for the n = 1 system this happens very quickly (cf. Table 
1^). There is a simple explanation of this phenomenon. As is noticed from the e.o.m. 
the iso-rotation with respect to the 1- and 2-axes 

fn{x) = 1 , 5'„(x) = c (zero-energy solution for a 0) (16) 

is a zero-energy solution of the M = n partial-wave for a — 0. The potential 
connected with a finite a ^ distorts this continuum solution into a bound-state, 
which in the limit a — is shifted towards the threshold. 



Table 2: Bound-state energies relative to the threshold in 
the M = n partial-wave. With decreasing parameter a the 
bound-states move towards threshold. 



1.0 



0.707 0.316 0.224 0.1 0.071 



n 

n 



0.943 0.99994 



1.0 



1.0 



1.0 



1.0 



0.330 0.394 0.622 0.690 0.866 0.943 



In the n = 1 system this is the only bound-state apart from the zero-modes 
already discussed. However in the n = 2 case some of the low-lying resonances 
may also become bound for sufficiently strong a, e.g. for a = 1 we find additional 
bound-states in the M = and M = 1 partial-waves at energies Uq = 0.923 and 
Ui = 0.793 respectively. All these bound-states must be considered for the Casimir 
energy (compare also the discussion of the phase-shifts in subsection 3.4). 
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3.3 Stability of hedgehog sohtons 

Formally, eqs. (|l^) are the e.o.m. for small amplitude fluctuations around a hedge- 
hog soliton with arbitrary topological quantum number n. For this reason they 
contain information about the stability of the n-skyrmion. For the ra-skyrmion to 
be stable there must not exist solutions to (13) with negative uj"^. 

Indeed, we do flnd solutions of (p!3| ) with negative uj"^ for topological quantum 
numbers > 3. For example for = 3 with = 0.1 there exist such solutions in 
the M = 2 (u;2 = -0.14), M = 3 (cu^ = -0.30) and aU higher partial waves. In the 
n = 1 and n = 2 sectors there exist no negative u"^ solutions. 

We conclude, that the n = 1 and n = 2 hedgehog solitons are stable while 
those with higher topological charges are unstable. This agrees with the flndings in 
ref. [|l^. In what follows we consider the scattering off the stable skyrmions with 
n = 1 and n = 2. 

3.4 Phase-shifts 



For the 2-dimensional scattering problem in spherical coordinates the incident 



plane wave is decomposed into partial waves carrying magnetic quantum numbers 

oo oo 

^ipr^ J2 t^^Jj,,{pr)e'^^=J2eMt''JMipr)cosiMip). (17) 

M=-oo M=0 

The latter transformation with the multiplicities eo = 1 and eM = 2 for M > 1 
follows from the properties of the regular and irregular Bessel functions of the flrst 
kind, Juiw) and NM^pr), when M — >• —M. Thus, it suflices to consider the partial 
waves M > 0. Similarly, for the scattered wave (plane wave plus outgoing radial 
wave) in the case of a single channel we have 

oo 

^ = ^Mi^'i'M{p,r)cos{Mif). (18) 

A/=0 

The partial-wave projected scattering waves iPm{p,^) with appropriate boundary 
conditions at the origin are integrated according to the e.o.m. and the phase-shifts 
then are obtained from the asymptotic form 

V'm(p, r) [JM{pr) cos 5m(p) + NM{pr) sin 5Mip)]e'''''^''^ . (19) 

The phase-shifts are related to the cross-section in the usual way, and the prop- 
erty Sm{p) = S~m{p) follows from the corresponding symmetry of the e.o.m.. The 
generalization to two coupled channels now is straightforward. It is noticed that 
our scattering equations (0) decouple asymptotically when the functions f± = 
(/m =F 9m) I ^ are introduced 

-;K)' + ^^^^/±=fV±, + (20) 
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with the solution 



f±{p, r) = A±{p)J\M±n\{pr) + B±{p)N\M±n\{pr) . (21) 

From the coefficients A±{p) and B±{p) of the asymptotical solution, the 2x2 scat- 
tering matrix Sm is obtained for every partial wave M. Again, as the e.o.m. are 
invariant with respect to M — >■ —M so is the scattering matrix and it suffices to 
consider M > 0. It is convenient to diagonalize the scattering matrix 

pi^M n \ 

Sm = Um\ \ ^,,.) (f/M)^ (22) 
e M I 

in order to obtain the intrinsic eigen phase-shifts 6^1^ and 6^^^ . Their sum 6m = 
"^M + is plotted for the partial waves with M < 4 in Fig. 3, where the n = 1 
and n = 2 systems are considered separately. 

The picture that emerges looks qualitatively quite similar to what has been 
obtained for 3D-SU{N) skyrmions long ago ||T3. For topological charge = 1 we 



obtain in detail: The phases for M = and M = 1, where the rotational and 
translational zero modes are located together with the M = 1 bound-state, start at 
TT resp. 2tt according to Levinson's theorem. For the smaller parameters = 0.01 
and 0.1 the wave-function of the M = 1 bound-state close to threshold (Table |]) is 
already quite similar to the infinitesimal translation (cf. subsection 3.5 for the limit 
a — >■ 0). As a consequence the phase-shift at threshold reacts with a sudden drop 
from 27r to vr which in the figures is only noticed because the phase-shift falls below 
TT before it bends to the right. A weakly pronounced breathing mode is observed 
at low energies for M = 0. The partial-waves M = 2, 3, . . . then contain a band of 
resonances which are flattened out for higher M and also with decreasing parameter 
a. 

Similarly for topological charge n = 2: For the smaller parameters = 0.01 and 
0.1 the M = and M = 1 phase-shifts start at tt because of the zero modes. In 
all but the M = 2 partial-wave with its bound-state (cf. Table ^ sharp resonances 
occur, which are much more pronounced compared to the n = 1 case. For strong 
potentials then, the M = and M = 1 resonances become bound (e.g. for a = 1 at 
ouq = 0.923 and uji = 0.793). Also less pronounced secondary resonances do appear. 



3.5 Belavin-Polyakov soliton 

Finally we would like to add a few comments on the pure Belavin-Polyakov solu- 
tion [jl2|. As mentioned we recover this solution for small parameters a with its 
size R still fixed by the balance of the Skyrme and potential terms (cf. appendix). 
In the pure Nia model without additional stabilizing terms the size R of the BP 
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topological charge n-1 
a=0.1 







topological charge n=2 
a=0.1 


















M=3 / M=4 /' 







topological charge n=1 






a=0.316 










M=1 








M=2 / ' 






M=3 













topological charge n=2 
a=0.316 




topological charge n=1 



topological charge n=2 
a=1 


M=0 


''m=i 




/ M=2 






/M=3 




M=4 





Figure 3: Phase shifts for the lowest partial waves scattered off the stable n = 1 
and n — 2 skyrmions. The relevant parameter was fixed at c? = 0.01, 0.1 and 1.0. 
The momentum p is plotted in natural units. According to Levinson's theorem the 
phase-shifts start at multiples of tt. Many resonances are observed in the various 
channels. 
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soliton is an undetermined parameter. In contrast to the full lagrangian (|1]) the pure 
Nia model does not break 0(3)-symmetry and is also conformally invariant. Due 
to these additional symmetries a larger number of zero-modes is expected. 

The e.o.m. for the fluctuations obtained in that case from (0) with a = 
and X = r/R decouple for f± = (/m =f gM)/V^ (not only asymptotically) and are 
easily solved. Still the counting of the zero-modes is a bit tricky. In addition to 
the rotation ([l^) , the breathing-mode (conformal invariance) becomes a zero-mode 
in the M = partial-wave as expected. But for the 1-soliton the iso-rotation 
(|16|) coincides with the translation (|T5|) giving rise to only 2x1 zero-modes in the 
M = 1 partial-wave. This makes altogether 4 zero-modes for the 1-soliton. Similar 
counting for the 2-soliton yields 8 zero-modes. 

Of course, this discussion is somewhat academic, because in reality there should 
always be some stabilizing term which fixes the size of the soliton and which at 
least breaks conformal invariance. Therefore we do not show the corresponding 
phase-shifts. 



4 Casimir energy 

With the bound-state energies and the phase-shifts provided in the previous section 
we are in the position to evaluate the 1-loop contribution to the soliton energy, i.e. 
the Casimir energy. The UV singularities contained in the loop are related to the 
high momentum behaviour of the phase-shifts 



M=0 



p >oo 2 

— y aop + ai . 



(23) 



In the case of a vanishing Skyrme term the coefficients 



ao = , 



ai 



+ 2m^(l -c) 



(24) 



are known analytically from the Born terms. For the full model they have to be 
extracted numerically from the phase-shift sum (^), which for that purpose has to 
be calculated to a high precision with some 100 partial waves taken into account. 
Typical values obtained, e.g. for = 0.1, are Oq = 1.0(/em)~^ and ai = 7.2 for 
the 1-soliton, and Oq = 2.0(/em)~^ and ai = 13.4 for the 2-soliton respectively. 
According to (^) the Nia plus potential term contributions to these coefficients are 
ai = 8.3 for n = 1 and ai = 15.6 for n = 2. With this established we may use the 
phase-shift formula f^, subtract the troublesome high momentum behaviour from 
the phase-shifts, and separately add the corresponding counterterms 



j^l—loop 

2n 



P 

dp—6{p) — m6{0) 

LU 



^ B 



(25) 
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This procedure is closely related to what has been employed in the 3 + 1 dimensional 



Skyrme model for the calculation of the Casimir energy to the nucleon mass pO| , |21 



The counter terms could e.g. be evaluated using a 2-momentum cutoff A >> m 



[A-m], / -^i- = — [A^--m'A + 2m^], (26) 



(27r)2 00 27v' " J (27r)2 oo Qtt' 2 

which makes the linear and cubic divergencies explicit 



j^l-loop _ ^0 

Qtt 



A' - -m^A 
2 



—A + E^"' . (27) 
2tt 



While the N^cx and the potential term contributions located in oi may be absorbed 
in a renormalized constant 

(/ denotes the bare constant), the divergencies stemming from the Skyrme term 
renormalize the Skyrme parameter e as well as the couplings of all the higher gradient 
terms not listed in the lagrangian (|l]). The latter renormalization we do not carry 
out explicitely, instead we simply assume that the renormalized Skyrme parameter 
be e and the renormalized couplings of the higher gradient terms be zero. To which 
extent this is a consistent assumption will be discussed in the following section 
in connection with the scale-dependence of our results. Now all the divergencies 
residing in ( P7| ) are taken care of, leaving the finite Casimir energy 

-yJ: < (*(^) - - - i - - 1^ . ^ . («) 

This result is independent of the employed regularization scheme. For instance, us- 
ing dimensional regularization the counter terms ( ^6]) in 2 + 1 dimensions become 
finite and equal to the last terms in the brackets, respectively, which yields the same 
expression for the Casimir energy. However, although dimensional regularization in 



odd dimensions has been used right from the beginning ^2[ and is now commonly 



applied in <l>3, QCD3 and Chern-Simons theories |^ we hesitate to apply this 



scheme for an odd number of dimensions because the fate of the UV singularities re- 
mains obscure. As the divergencies do not show up in odd dimensional regularization 
there seems to be no possibility to introduce a renormalization scale. For strictly 
(super) renormalizable theories this may be of no further importance, however in 
theories which are renormalized order by order in terms of a gradient expansion a 
scale is desirable in order to control the convergence of the series (cf . section 5) . 
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-4 I 1 1 ' 1 

12 3 4 

momentum p 

Figure 4: Subtracted phase-shift sum as a function of the momentum p in natural 
units for the n = 1 and n = 2 systems calculated with = 0.1. 

The Casimir energy (^91) consists of three parts, i.e. the phase-shift integral, the 
bound-state contribution and the contribution from the counter-terms. 

The subtracted phase-shift sum S{p) — a^p^ — ai which enters the phase-shift 
integral is plotted in Fig. 4 for the n = 1 and n = 2 systems {a? = 0.1). The 
maximum momentum to which this sum has to be integrated depends sensitively 
on the parameter a. The various contributions to the Casimir energy are given 
in Table ^ for several values of a. The dependence of the Casimir energies on 



Table 3: Individual contributions to the Casimir energy according 
to ( p9D for the n = 1 and n = 2 systems. All energies are given in 
natural units. 





n = 1 


n = 2 


a 


0.1 0.316 1.0 


0.1 0.316 1.0 


phase-shifts 
bound-states 
counter terms 


-0.29 -0.08 0.00 
-0.47 -0.84 -1.56 
0.31 0.62 1.85 


-0.63 -0.26 0.05 
0.52 -1.06 -2.47 
0.61 1.16 3.36 


total 


-0.45 -0.30 0.29 


-0.54 -0.15 1.00 



this parameter is shown in Fig. 5. Note that for an easier comparison again only 
half of the 2-skyrmion's Casimir energy is plotted. It is noticed that the Casimir 
energy of the 1-skyrmion stays always smaller compared to half of that of the 2- 
skyrmion. Therefore, in contrast to the classical energy the Casimir energy favors 
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Figure 5: Casimir energies in natural units for the 1- and 2-skyrmion devided by 
their topological charge. 



single skyrmions. These competing effects may be summarized in the formula for 
the total energy cast into the form 



E, 



total 



AnfE'M + ^femElia) 

ATTf lE'M+yKia) 



y 



\/ fern 
1^ 



(30) 



where E\{a) represents the Casimir energy (pOf) in natural units \fjem. Thus, apart 
from the overall scale 47r/^ the model is characterized by the two dimensionless 
parameters a and y. If this ratio y exceeds a critical value (which increases slowly 
with a, e.g. y = 0.3,0.4,0.7 for = 0.01,0.1,1.0), then Ef*"' - 2Ef*"' becomes 
positive and the 2-skyrmion decays into two single 1-skyrmions. 



4.1 Numerical example from hadron physics 

Although different dimensionalities may cause qualitative differences, let us illustrate 
these results by a numerical example using the scale of hadron physics: Take the 
value a = 0.1, fix the size scale at Xj^fjem = 0.2 fm (y/ fern = 1 GeV) and make 
use of the numbers given in Table 1 and Table 3. With the above parameters the 
1-soliton gets a topological radius 0.27 fm and a Casimir energy E^'^'^ = — 0.45GeV, a 



situation very similar to what has been obtained in the hadronic Skyrme model |21 
if in addition the classical soliton mass is fixed at Ef"'^^ = 1.21(4??/^) = 1.39 GeV 
in order to obtain the nucleon mass at El"^"^ = 0.94 GeV. Then the n = 2 soliton's 
total energy E^"^""^ = (2.64 - 0.53) GeV = 2.11 GeV turns out to be larger than 
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that of two separated 1-skyrmions. The parameter in this example, y = 0.9, hes far 
above the critical value 0.3. Of course, the 2 + 1 dimensional model discussed here 
cannot be taken literally for hadron physics, but exactly this mechanism may shift 
the undesired torus configuration obtained in hadronic soliton models with baryon 
number B = 2 to higher energies. 



4.2 Casimir energy of the pure Belavin-Polyakov soliton 

Finally we would like to discuss the Casimir energy of the pure BP soliton with free 
size parameter R in order to make contact to the quantum corrections as obtained 
in ref. |^ . With m = and the coefficients = , ai = 27m and no additional 
bound-states present (cf. the discussion about the zero-modes in subsection 3.5) 
eq.(|25|) for the 1-loop contribution simplifies 

(A) 



j^^.loop BP _± n >^^^^^^ 
271 Jo 



271 Jo 

^,^'^^dp{5{p)-a,)-^K (31) 

ZTT JO 271 

1 rim 



/ d{pR) {6{pR) - 27m) - nA . 
Jo 



2nR Jo 

We indicate here the possibility to limit the momentum integration by a finite Debye 
momentum A = p£, related to the lattice constant. Because the soliton should at 
least cover several lattice sites, AR > 1 is large enough to extend the integration over 
the subtracted phase-shifts to infinity (see Fig. 4) just as in our field theoretical 
treatment. Through the substitution in the last step the dependence of the sub- 
tracted phase-shift integral on the soliton size R becomes explicit: the magnitude 
of the Casimir energy decreases with increasing soliton size. 
Numerically, we find 

for the pure BP solitons with topological charges n = 1 and n = 2 respectively. Their 
sizes Ri and R2 are arbitrary and generally different. If the sizes are determined by 
the stabilizing terms in the limit a ^ 0, we obtain _R„ = 6„/ \fjem with 61 — > and 
62 — ^ "\/2 (appendix). For that reason the Casimir energy of the 1-soliton plotted 
in Fig. 5 bends to negative values and finally diverges as a approaches zero. The 
Casimir energy of the 2-soliton stays finite, E\ = —1.7/\/2 + 0{^/a), but with an 
infinite slope at a = 0. 

In ref. the subtracted phase-shift integral was neglected for pdR >> 1 and the 
last term in ( PTf ) then lead to the result E^''^°°p ^ —npi). Because this is exactly the 



term we absorb into a renormalized / (p7|, 28 ) , this result corresponds to E'™* = 
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in our notation. The opposite limit pnR << 1 where the 1-loop contribution (|3T| ) 
vanishes corresponds to sohton sizes much smaller than the lattice spacing and does 
not make sense. The conclusion, that the quantum corrections lower the soliton 
energy as soliton size increases, is wrong. 



5 Scale dependence 

In this section we introduce a scale /i which allows to shift finite pieces from the tree 
to the 1-loop contribution and vice versa. It will be introduced in such a way that 
we recover the results of the previous section for fi = fiQ. Tentatively, /io = 47r/^ 
may be identified with the underlying energy scale of the model in analogy to the 
chiral scale /io = 47r/7r ~ IGeV of ChPT ||2^. Our results will however be presented 
in such a way that no fixation of fiQ is required. For this purpose we write the 
cutoff characterizing the divergencies of the 1-loop contribution (^) as a sum of 
two parts A = [A — (yU — yUg)] + (/i — fio)- The first part is then renormalized into a 
scale-dependent strength of the N£cr and potential term 

f (/i) = f - ^-^^-^'^ =f + . (33) 

According to this formula fi ^ fio shifts the energy scale 47r/^(/i) away from its 
original value 47r/^. In order to exhibit the scale-dependence it is now convenient 
to present our results as function of the ratio /^(/u)//^ rather than /i itself. The 
presentation is then independent of /iq and also robust against changes in the regu- 
larisation scheme which may introduce numerical factors to the scale (e.g. in the 
3-momentum instead of the 2-momentum cutoff scheme the scale is changed by a 
factor of 7r/2). The second part of the above decomposition, (/i — /xq), remains 
explicit in the finite Casimir energy 



E™^(/x) 



/TT JO UJ 



aop^ - ai) + ^ ^{ujb - m) 

^ B 



(34) 



+ 



Oo 

Gtt 



(// - /io)^ 



3m^ 



(/X - /xo) - 2m^ 



2^ 



[/X - /xq - m] . 



Also for the Qq - term with its cubic divergence we kept all contributions from the 
second part in the Casimir energy although other prescriptions are possible. However 
this ambiguity will not infiuence the results strongly because of the smallness of 
Qq (cf. numbers in the previous section). In this way both, the tree and the 1- 
loop contributions become scale-dependent. While for the N£o" and potential terms 
the scale-dependence in the total soliton energy is exactly compensated, the 1- 
loop contribution which arises from the Skyrme term would not only yield a scale- 
dependent Skyrme parameter e(/x), but would also switch on all higher gradient 
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terms assumed to be zero for jj, = hq. Having assumed all these couplings to be 
independent of the scale implies that the tree + 1-loop contribution to the soliton 
energy cannot be strictly scale-invariant. In fact, the resulting scale-dependence 
measures the magnitude of the higher gradient terms not accounted for. All the 



more it comes as a surprise, just as in the hadronic case that for specific 
parameter choices an almost scale-independent soliton mass may be obtained. Such 
a case (a^ = 0.1,?/ = 1) is depicted in Fig. 6 for the 1-soliton. The rather strong 
scale-dependence of the tree contribution is nicely compensated over a wide range 
P{lj) = 0.5p, . . . ,2.0p when the 1-loop contribution is added. Of course, this 
statement has to depend on the parameters used. Therefore in Fig. 7 we plotted 
the total 1-soliton mass for = 0.1 and various values of the ratio y. It is noticed 
that weak scale-dependence requires a ratio close to ?/ ~ 1, the value used in the 
previous figure. Lower and higher values of y enhance the scale-dependence. We do 
not show the corresponding plots for the 2-soliton because they look quite similar, 
with the scale-dependence being even somewhat weaker in that case. 

We may now pose the question for which parameter combinations a and y we 
may expect a modest scale-dependence in accordance with the assumption that the 
couplings of all higher gradient terms are zero. The answer is illustrated in Fig.8. 
Inside the inner contour the average scale-dependence in the intervall P{^Jt)/ P G 
[0.5, 2] is less than 3%. For comparison the 6% contour is also plotted. 

All the parameter combinations lying inside the contours yield negative Casimir 
energies whose absolute values though sizeable do not exceed 50% of the classi- 
cal soliton mass. They also lead to an unstable 2-soliton which decays into two 
individual 1-solitons. In particular the parameters (a = 0.1,?/ = 0.9) of our nu- 
merical example from hadron physics (section 4.1) lie in the center of the almost 
scale-independent region. 

The weak scale dependence obtained for parameter combinations lying inside the 
contours in Fig. 8 implies that loop-corrections can be reliably calculated within 
the framework of the lagrangian (|^) with all higher couplings set equal to zero. 



6 Conclusions 

We studied the magnon-vortex system in the 2-1-1 dimensional 0(3) model in a field 
theoretical approach. For that purpose the N^cr model was augmented by stabilizing 
standard fourth-order Skyrme and potential terms. 

Complete information on bound and scattering states in all partial-waves was 
established. We find an extremely rich excitation spectrum with pronounced res- 
onances as in the ?>D-SU{Nf) case with its baryon resonances. In principle these 
resonances should be accessible by measuring the excitation spectrum of spin-waves 
in the presence of skyrmions. 

Furthermore, the quantum fluctuations will allow to study their influence on the 
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Figure 6: Scale-dependence of the 1-soliton's energy in tree and tree + 1-loop. The 
dimensionless parameters are a — a(/io) = -316 and y = |/(/Xo) = 1.0. The energy is 
in units 47r/^ and the scale is given by the ratio P{lJ)/ P- 
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Figure 7: Scale-dependence of the 1-soliton's energy in tree + 1-loop for a — a{iJ,o) — 
.316 and various parameters y = y^iio). Energy and scale as in Fig. 6. 
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Figure 8: Parameter space (a, y) for 1-skyrmions. Inside the inner contour the 
scale-dependence is less than 3%, inside the outer contour less than 6%. 



shape of the soliton Note that in tree approximation the third component of the 
order parameter field is always tied to —1 at the origin for topological reasons, in 



variance with microscopic Hartree-Fock calculations for the ferromagnet p6[ which 
suggest a vanishing value in case of small solitons. 

Finally with the bound state energies and the scattering phase-shifts the Casimir 
energies were evaluated. The appearing UV singularities were renormalized under 
the assumption that the (renormalized) couplings of the higher gradient terms be 
small. A criterium for the consistence of this assumption is the approximate scale- 
independence of the results such that the effective action represents an almost 1-loop 
renormalizable theory. This criterium requires a detailed balance of tree and 1-loop 
contributions which limits the parameter space where this requirement is met. The 
following results apply for this restricted parameter space. Parameters lying outside 
that range may lead to other conclusions, but there the 1-loop contribution cannot 
be reliably calculated within the model. 

• The Casimir energy is negative and generally large leading to a considerable 
reduction of the total soliton energy. However, it does not exceed 50% of the 
tree contribution and thus may still be considered a (sizeable) correction. 

• The magnitude of the Casimir energy decreases with increasing soliton size. 
But it also decreases as the stabilizing terms become more important relative 
to the N£(j model term. 

• With the 1-loop corrections included the n = 2 soliton becomes unstable and 
decays into two individual n = 1 solitons (which may still be weakly bound 
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by a dipole force). We conjecture that the same situation might occur for the 
3D-SU{2) B = 2 torus in hadron physics. 

It is obvious that the program outhned in this paper has to be tailored for specific 
apphcations concerning ferromagnets and antiferromagnets. Most importantly for 
ferromagnets the time derivative part of the lagrangian has to be replaced by the 
T violating Landau-Lifshitz term with one time derivative only. This replacement 
may change the results for the Casimir energies appreciably. Also, for antiferromag- 
nets an external magnetic field should couple through the time component of the 
covariant derivative. Finally, proper consideration of the non-local Coulomb inter- 
action complicates the situation: instead of coupled differential equations coupled 
integro-differential equations have to be solved for the fluctuations. 

The evaluation of loop corrections as presented here naturally suggest the in- 
clusion of temperature into the formalism. This opens the possibility to study the 
properties of 2D spin-textures at finite temperature, which is of particular interest 
near the symmetry restoring phase-transition. 

This work is supported in parts by funds provided by the FCT, Portugal (Contract 
PRAXIS/4/4. l/BCC/2753) . 

Appendix 

In this appendix we present analytical soliton solutions for the two limiting cases 
a — s> and a — oo. 



Belavin-Polyakov solution for a ^ 

The analytical solution of the Euler-Lagrange equation (|^) in the limit a — 

-(a;F')'-^ = (A.l) 
with boundary conditions (|^) is the well-known Belavin-Polyakov soliton 

ta4^(M". (A.2) 



The size parameter 6^ = y2n(n2 — l)/3 is determined by the interplay of the 
Skyrme- and potential terms contributing to the static energy functional 



K = n + a — -. - + a^-^, n^l. (A.3) 

3o„ smir/n 2 sm it /n 
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The corresponding square radius follows from its definition (^) 



<x2>„=W— 722-1 —X^ riy^l. A.4 
V 3 smiT/n 

It is noticed that in case of the 1-soliton the contribution from the potential term 
to (A. 3) diverges. For that reason we modify (A. 2), 



F / e 

tan — = 6.n W — 5 , (A. 5) 

2 V e^^ - 1 ^ ^ 

such that for e — the Belavin-Polyakov solution is recovered. The variation of 
the static energy functional in lowest order e 

gives then e = and ine = —A/Sbf with the result 

^0 = 1 + e(i - £ne) , = —/ine . (A.8) 

Indeed with a also e tends to zero. The square radius behaves like 



< >i= y -^^"-e (A.9) 

and diverges weakly in the limit a — 0. For a ^ 0.001 the radius of the 1-soliton 
exeeds that of the 2-soliton. Exactly this behaviour may be traced in the numerical 
solution if the differential equation (|^) is solved for very small a with great precision. 
Because the square radius (A.9) equals the derivative of the soliton energy with 
respect to the parameter a, this also explains the infinite slope at the origin indicated 
in the 1-soliton's curve plotted in Fig. 2. 

Solution for a — oo 

The Euler-Lagrange equation (|^) in the opposite limit a — » oo 

X \ X ) 

may also be solved analytically 



-s = (A.IO) 



cosF(x) = -^(62 -x2y + 2^-l, x<6„. (A.ll) 
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The integration constants are fixed by the condition F{bn) = connected with 
topological charge n. The size parameter 6^ = 4n again minimizes the static energy 
functional 



= na 

with the results quoted in (|^) 



An 2\bl) 6\An) 



(A.12) 



E0 = — a, <a;2>„=_. (A.13) 

Needless to say, that numerical solutions of the Euler-Lagrange equation (0) repro- 
duce these results for large enough a. 
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